Ground states of trapped spin-1 condensates in magnetic field 
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We consider a spin-1 Bose-Einstein condensate trapped in a harmonic potential under the influence 
of a homogeneous magnetic field. We investigate spatial and spin structure of the mean-field ground 
states under constraints on the number of atoms and the total magnetization. We show that the 
trapping potential can make the antiferromagnetic condensate separate into three, and ferromagnetic 
condensate into two distinct phases. In the ferromagnetic case, the magnetization is located in the 
center of the harmonic trap, while in the antiferromagnetic case magnetized phases appear in the 
outer regions. We describe how the transition from the Thomas-Fermi regime to the single-mode 
approximation regime with decreasing number of atoms results in the disappearance of the domains. 
We suggest that the ground states can be created in experiment by adiabatically changing the 
magnetic field strength. 

PACS numbers: 03.75.Mn, 03.75.Hh, 67.85.Bc, 67.85. Fg 
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I. INTRODUCTION 

Bose-Einstein condensates with spin degrees of free- 
dom [lj attracted in recent years great interest due 
to the unique possibility of exploring fundamental con- 
cepts of quantum mechanics in a remarkably controllable 
and tunable environment. The ability to generate spin 
squeezing and entanglement [2] makes spinor Bose gases 
promising candidates for applications as quantum sim- 
ulators [3|, in quantum information [J], and for precise 
measurements [5J- Moreover, spinor condensates were 
successfully used to recreate many of the phenomena of 
condensed matter physics in experiments displaying an 
unprecedented level of control over the quantum system. 
In particular spin domains [H-H, spin mixing and 
spin vortices [llj were predicted and observed. 

The ground states of spin-1 condensates in homoge- 
neous magnetic field have been studied in a number of 
previous works [f| Ii~2r{l6| . The typical procedure consists 
of minimization of the total energy under constraints on 
the number of atoms and the longitudinal magnetization, 
which is well conserved in typical experimental condi- 
tions Most of the studies, however, considered the 
condensate in the single-mode approximation, which as- 
sumes that the spin components share the same spatial 
profile [H], 03 j ignoring the possibility of phase separa- 
tion. 

In Ref. [IH , the breakdown of the single-mode approx- 
imation was shown numerically for a condensate confined 
in a harmonic potential. In a recent paper [l6| . we stud- 
ied the possibility of phase separation in the case of an 
untrapped condensate, and demonstrated that this phe- 
nomenon can take place in an antiferromagnetic conden- 
sate, leading to the splitting into two distinct phases. On 
the other hand, ferromagnetic condensates were found 
not to display phase separation for any values of pa- 
rameters. However, the trapping potential, present in 
any experimental realization, can strongly influence the 
structure of the ground states. Indeed, in the context of 
binary condensates, it was previously shown that spatial 
separation of components can be induced by an external 



potential This phenomenon was called "potential 

separation" as opposed to "phase separation" which oc- 
curs spontaneously also in the untrapped case. Moreover, 
the external potential determines the spatial structure of 
phase separated states. In general, phases that are char- 
acterized by energetically unfavorable interatomic inter- 
actions are moved towards the outer, low density regions, 
as demonstrated in spin-imbalanced Fermi systems [l9j 
and spinful bosons [Til H3| ■ 

In this paper, we investigate in detail the spatial sepa- 
ration of spin phases in the ground states of spin-1 con- 
densates trapped in harmonic potentials. We find that 
the phase diagram of a trapped condensate is substan- 
tially different from the one of an untrapped condensate. 
We show that the external potential can make the anti- 
ferromagnetic condensate separate into three, and ferro- 
magnetic condensate into two distinct phases. Further- 
more, we suggest an experimental method for the cre- 
ation of these ground states by adiabatically changing 
the magnetic field strength. 

The paper is organized as follows. Section |TT] reviews 
the mean-field model of a spin-1 condensate in a ho- 
mogeneous magnetic field and the possible spin phases. 
Section ITVl presents analytical and qualitative results de- 
scribing the structure of ground states in the Thomas- 
Fermi approximation. In Section [V] we treat the problem 
numerically in a systematic way and demonstrate the 
crossover between the Thomas-Fermi and single-modc- 
approximation regimes. Section IVII describes an exper- 
imental method of creation of the ground states. Sec- 
tion |VII] concludes the paper. 



II. MODEL 



We consider a dilute spin-1 BEC in a homogeneous 
magnetic field pointing along the z axis. We start with 
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the mean-field Hamiltonian H 

h 2 _„ 



H S + H A , 



H 



3=-,0,+ 



2m 



Co 
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n + V(r)^ ^+F 4 , 



where ?/>_ , tpo > "0+ are the wavefunctions of atoms in mag- 
netic sublevels rrif = — 1,0, +1, M is the atomic mass, 
V(r) is an external potential and n = J2 n j IV^'I 2 is 
the total atom density. The asymmetric (spin dependent) 
part of the Hamiltonian is given by 



H A = / drl ]T 



j=-,o,+ 




(2) 



where Ej is the Zeeman energy shift for state tpj and the 
spin density is, 

F = (F x ,F y ,F z ) = (^F^^Fy^^F^) (3) 

where F XtVtZ are the spin-1 matrices [21] and ip = 
(ip+, "00, V'-)- The spin-independent and spin-dependent 
interaction coefficients are given by cq = ATih 2 (2a 2 + 
ao)/3m and c 2 = Airh 2 (a 2 — ao)/3m, where a$ is the 
s-wave scattering length for colliding atoms with total 
spin S. The total number of atoms and the total magne- 
tization in the direction of the magnetic field 



N 



ndr . 



M = J F z dr = J (n+ — n_) dr , 



(4) 
(5) 



are conserved quantities. The Zeeman energy shift for 
each of the components, Ej can be calculated using the 
Breit-Rabi formula (23| 



E± = --£ H FS (l + 4 V / l±a + a 2 ) qp 9iPbB , 



En = — — E 



HFS 



1 + 4^1 + a 2 



(6) 



where -Ehfs is the hyperfine energy splitting at zero 
magnetic field, a — {gi + g.])pBB / E^ys, where pb is 
the Bohr magneton, gj and gj are the gyromagnetic ra- 
tios of nucleus and electron, and B is the magnetic field 
strength. The linear part of the Zeeman effect gives rise 
to an overall shift of the energy, and so we can remove it 
with the transformation 

H -> H + (N + M)E + /2 + (N - M)E_/2. (7) 

This transformation is equivalent to the removal of the 
Larmor precession of the spin vector around the z axis 
fl~3l [23j . We thus consider only the effects of the 
quadratic Zeeman shift. For sufficiently weak magnetic 
field we can approximate it by SE = (E + +E- —2Eq)/2 
a 2 £ , HFs/16, which is always positive. 

The asymmetric part of the Hamiltonian ^ can now 
be rewritten as 



// a = / <lr M£n + y|F| 2 ) = 



drne(r), (8) 



where the energy per atom e(r) is given by pj 
e = -5E PQ + ^f\i\ 2 = -SE P0 + C -f(\{ ± 



\ 2 +m 2 ) 



(1) |f ± | 2 =2po(l - Pa) + 2^0^(1 - Po) 2 - m 2 cos(0) . (9) 



We express the wavefunctions as ijjj 



/npj exp(i9j) 



where the relative densities are pj = rij/n. We also in- 
troduced the relative phase 6 = 9 + + 6- — 29n, spin per 
atom f = F/n, and magnetization per atom m = f z = 
p + — p- . The perpendicular spin component per atom is 

\^\ 2 = fl + fl 

The Hamiltonian ([1} gives rise to the Gross-Pitaevskii 
equations describing the mean-field dynamics of the sys- 
tem 

ih ~tt~ = ^ + c z( n ± +n - n T )} ip± + c 2 ip 2 ip^ , (10) 
dip 

ik-Qj- = [C - 5E + c 2 (n + + n-)} ip + 2c 2 -)jj + ip-')pQ , 

where C is given by C — —h 2 V 2 /2M + cnn + V(r). 

By comparing the kinetic energy with the inter- 
action energy, we can determine the healing length 
£ = 2irh/ yJ2Mc§n and the spin healing length £ s = 
2irh/ 1 \J2Mc 2 n. These quantities give the length scales 
of spatial variations in the condensate profile induced by 
the spin- independent or spin-dependent interactions, re- 
spectively. Analogously, we define the magnetic healing 
length as £ B = 2Trh/v / 2M5E. 

In spin-1 condensates created to date, the ao and a 2 
scattering lengths have similar magnitudes. The spin- 
dependent interaction coefficient c 2 is then much smaller 
than its spin- independent counterpart cq. For example, 
this ratio is about 1:30 in a 23 Na condensate and 1:220 
in a 87 Rb condensate far from Feshbach resonances [2~jj . 
Consequently, changing the total density n requires much 
more energy than changing the relative populations of 
spin states rij. In our considerations we will treat the 
total atom density profile n(r) as a constant, close to the 
Thomas- Fermi profile for a given potential V(r). 



III. HOMOGENEOUS STATIONARY STATES 

We recall the possible phases of spin-1 condensates in 
magnetic field [l6| . The stationary solutions in the case 
of a vanishing potential, V(r) = 0, have the form 



ipj(r,t) = ^ 



(11) 



where ps — cqTi/H is a constant. These solutions are sta- 
tionary in the sense that the number of atoms in each 
magnetic sublevel rij is constant in time, but the rela- 
tive phases may change as a result of an additional spin 
precession around long as the phase matching con- 
dition 



p- = 2p , 



(12) 
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is fulfilled [23J, l25|. Because the symmetric part of the 
Hamiltonian in ([1]) is constant, the relevant part of the 
Hamiltonian is given by Eq. (jHJ. 

The Hamiltonian (p} and GP equations (fTUf are invari- 
ant under the gauge transformation ipj — > ipj exp(— i/3) 
and rotation around the z axis ipj — > ipj exp(—iF z ~/), 
which transform the wavefunction components according 
to 



(13) 



Hence the solutions can be classified using the rela- 
tive densities pj — n,j/n and a single relative phase 
9 = 9 + + 6- — 26*o, with the chemical potentials fij given 
as solutions to Eqs. (fTU)) . We note that for stationary so- 
lutions with all three components populated the relative 
phase must take one of two values, 9 = or 9 = ir. We 
call the former "phase-matched" states (PM) and the lat- 
ter "anti-phase-matched" (APM) states. The names de- 
rive from the fact that within the continuum of states sat- 
isfying the spin rotations (| 13[) there is a set (ip+,ipo,ip-) 
with all components in phase for the PM states, and a 
set with ip + and ipo in phase but with ip_ ir out of phase 
for the APM states. 

In general, the possible stationary states can be clas- 
sified as follows 

1. Nematic state (po), with all the atoms in the to/ = 
component, tp — (Q,ipo,0). The chemical poten- 
tial, spin density and energy per atom are equal 
to 



hp = -5E, |f|=0, e=-SE. 



(14) 



2. Magnetized states (/?_ and p+), with all the atoms 
in the m,/ = —1 or in the to/ = +1 component, 
^ = (0,0,V+) or^ = (^_,0,0) 



hp± = c 2 n, \m\ = 1, |f]_| = Q, e =2' ( 15 ^ 



3. Two-component or stretched states (2C), with ip = 
(ip-, 0, ip+) and arbitrary magnetization 

hp± = ±C2?ito, |fj_| = 0, e — • (1^) 

4. Phase-matched states (PM), where all three sub- 
levels are populated, ip — (V'-j rfoi ip+)> an d the 
relative phase is equal to 6 = 0. 

5. Anti-phase-matched states (APM), similar as the 
PM state but with 9 = ir. 

The parameters of the PM and APM phases cannot 
be expressed by analytical formulas in the general case. 
Since all three components are populated in these states, 
the magnetization perpendicular to the magnetic field 



direction, |Fj_| = v / 2|V , + 1 /'o + V'oV'-li is nonzero. Hence 
the axial symmetry of the system is broken [l3[ . 

These above classification is also applicable to inho- 
mogeneous condensates within the single-mode approxi- 
mation (SMA), which assumes that the spin components 
share the same spatial profile [3, [ljj , after replacing n 
with (n). This assumption is true eg. when the conden- 
sate size is much smaller than the spin healing length £ s 
and the magnetic healing length £b. 



IV. GROUND STATES IN THE 
THOMAS-FERMI APPROXIMATION 

In this section, we investigate the structure of the 
ground states in harmonic potential within the Thomas- 
Fermi (TF) approximation (or local density approxi- 
mation), neglecting the kinetic energy in the Hamilto- 
nian |T]). This approximation is justified when the spatial 
variation of the condensate wavefunction gives a contri- 
bution that is relatively small, in particular when the size 
of the condensate is much larger than both ^ s and 

Within our assumptions, the asymmetric part of the 
energy © is a small contribution to the total energy. 
The profile of the total density is determined by the min- 
imization of the symmetric part of the free energy with 
a fixed particle number 



V )n dr - p,N , 



(17) 



where p, is the Lagrange multiplier corresponding to the 
chemical potential. The minimization of this functional 
gives the well known Thomas-Fermi profile con(r) + 
V(r) = p, = const. 

Analogously, to determine the spin structure, we min- 
imize the asymmetric part of the energy Ea with a fixed 
total magnetization M 

Ka = J e(n, m)n dr — pM , M = J mndr. (18) 



The minimum of this functional requires that 
de(n, to) 



dm 



p = const, 



(19) 



when moving from one point in space to another. This 
condition is not applicable to the completely magnetized 
p+, p_ and the unmagnetized po phase, for which de/dm 
is undefined 11611 - 



A. Antiferromagnetic case (c2 > 0) 

Since we work within the Thomas-Fermi approxima- 
tion and neglect the spatial derivatives, at each point in 
space the local phase corresponds to a ground state of a 
homogeneous condensate [la |. Nevertheless, the change 
of local parameters (density and magnetization) in space 
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FIG. 1: Schematic profiles of phase separated states, (a) 
The swapping of two domains can change the total energy, 
even if the number of atoms in each of them is kept constant, 
(b) The structure of the ground state of an antiferromagnetic 
condensate composed of three phases. The dashed line shows 
the local magnetization. 



may lead to the change of the local ground state and the 
appearance of spatial domains. 

In the case of an antiferromagnetic condensate, there 
are five possible local ground state phases: p_|_, p_, po, 2C 
and APM [l(| . The last two of them allow for a variation 
of the local magnetization. In the case of the 2C phase, 
the condition (|19l) gives 



de 
dm 



const, 



(20) 



hence the difference in the density of atoms in the m = 
+1 and the m = — 1 components is constant. For the 
APM state, the asymmetric energy (further on we call 
it simply the energy) is not an analytical function of n 
and m, and we approximate it with the formula eAPM ~ 
(|m| — 1)6E+ \m\c2n/2, which gives less than 10% error 
for any n and m. The condition (|19[) gives 



de 
dm 



5E 



c 2 n 



const, 



(21) 



which obviously cannot be fulfilled if the density profile 
is not homogeneous. Consequently, the APM phase is 
unstable, because the transfer of magnetization from the 
central, high density area to the outer regions is always 
energetically favorable. Hence, the ground state can con- 
tain p + , p_, pq, and/or 2C domains, the latter fulfilling 
the condition (|20p . 

The other question is the placement of the domains. 
Since the local ground state depends on the density, the 
domains arrange according to the distance from the cen- 
ter of the trap, as in Fig. QJa). Let us consider two do- 
mains, one of them composed of phase X placed at the 



area 1 and the other composed of phase Y placed at the 
area 2, containing equal number of atoms. If the energy 
per atom of the phase X is constant, (which is the case 
for the po phase), then swapping of the domains as in 
Fig. [TJa) will not lead to the change of its energy 



AE X = e x 



n dr 



n dr 



0. 



(22) 



On the other hand, if the energy per atom of the phase 
Y is proportional to the number of atoms, ey = an, as 
for the p± or 2C phases, the difference in its energy is 



n 2 dr — / n 2 dr) = aB, where B > . 



\Jl J2 1 

(23) 

In result, it is energetically favorable for the po phase to 
stay in the high density area in the center of the trap. 
Moreover, the p+ or p_ phase will be moved further than 
the 2C phase to the outer areas, since their coefficient a 
is larger. 

We obtain the general domain structure depicted in 
Fig. Hlb), shown for positive M (for negative M, the 
local magnetization is opposite and the p + state is re- 
placed by p- ) . The boundary between the phases po and 
2C corresponds to a first order phase transition, and the 
boundary between the phases 2C and p+ to a second or- 
der phase transition jig]. We note that not all of the 
phases must be present. If all three domains are present, 
we can derive the relation between n\ and n 2 by cal- 
culating the change of the total energy when changing 
the domain sizes a, b and c. Since the magnetization is 
conserved, we have AM = bAn 2 — n 2 Aa = and conse- 
quently An 2 /n 2 = Aa/b. The change of energy is 



AE A = (-6Eni - y^J Aa + yA(^) 



(-2 



-SEm 



C 2 2 



n Aa. 



(24) 



This gives the equilibrium condition 5En\ — c 2 n^/2, and 
consequently the necessary (but not sufficient) require- 
ment 



SE < yn max 



(25) 



for the coexistence of all three phases, where n max is the 
maximum density in the center of the trap. 



B. Ferromagnetic case (c2 < 0) 

In the case of a ferromagnetic condensate, the possible 
ground state phases are p+, p-, po, and PM 16]. Since 
it is difficult to propose a simple analytical formula for 
the energy of the PM phase, we resort to a qualitative 
analysis using the condition (fl9|) . 

In Fig. [2] we show the normalized energy per atom 
e/SE of these states in function of the local magneti- 
zation m for several values of the normalized density 




e/5E 

FIG. 2: Normalized energy per atom versus magnetization 
for the possible ground states of the ferromagnetic conden- 
sate. The lines correspond to normalized densities (from top 
to bottom) \c 2 \n/SE = 0.1, 0.4, 0.7, 1.0, 1.3. All the paths ful- 
filling the condition de/dm — const, with the exception of 
completely magnetized states, arrive at the central point cor- 
responding to the po state in the low n limit. 
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FIG. 3: Examples of ground states of an antiferromagnetic (a- 
c) and a ferromagnetic (d) condensate in the quasi-lD case, 
(a) 2C + p+ state for m — 0.8 and <5-E/c2n max = 0.05, (b) 
po + 2C + p+ state for m = 0.6 and <5-E/c2?i max = 0.17, (c) 
po + p+ state for m = 0.4 and 5E/c2n ul ^ — 0.37, and (d) 
PM + po state for m — 0.1 and J-E/c2?imax = —1.13. The n+, 
no, and n_ components are depicted by dash-dotted, dashed, 
and dotted lines, respectively, and the solid lines show the 
total density. Here u>± — 2tv x 10 3 Hz, other parameters are 

LUz = 2ir x 2.5 Hz in (a-c) and 
bj z = 2tt x 1.7Hz in (d). 



N = 8.4 x 10 J 



Na atoms 



N = 8.4 x lO 3 87 Rb atoms, 



c^ji/bE. The central point corresponds to the po phase, 
while the p_|_ and p_ phases correspond to the limits 
to = ±1. The possible Thomas-Fermi profiles correspond 
to sets of points in the e, to plane that fulfill the condition 
de/dm = const. At constant SE, the magnitude of this 
derivative generally increases with decreasing n (when 
moving from the center of the trap to the less populated 
regions), which can be seen from the gradually "steeper" 
shape of the top curves. To keep de/dm constant, the 
magnitude of magnetization has to decrease, and even- 
tually we arrive at the point po at a low enough n. In 
this way, we obtain a phase separated state consisting of 
PM and po domains. The only exceptions from this are 
the purely po state for M = and high enough SE, and 
purely p + and p_ states for M = ±1. 

The placement of domains in the PM+po state can 
be deduced in a similar way as in the antiferromagnetic 
case. The energy per atom in the po phase is constant, 
while the energy in the PM phase is a decreasing function 
of density, due to the spin-dependent atomic interaction 
with a negative C2 . Hence the atoms in the PM state will 
reside in the trap center, where the energy is the highest. 



V. NUMERICAL RESULTS 

A. Potential induced spatial separation 

In this section we present examples of numerically de- 
termined ground state profiles for both antiferromagnetic 
and ferromagnetic condensates in quasi ID and quasi 2D 
setups. The ground state profiles for a quasi-lD conden- 



sate were found numerically by solving the ID version 
of Eqs. (fTO]) [16l | with rescaled interaction coefficients 
c o D ' c 2 D = ( muj ±)/(^ 7T ^) c 0i c 2, where uj± is the trans- 
verse trapping frequency. The Fermi radius of the trans- 
verse trapping potential is smaller than the spin healing 
length, and the nonlinear energy scale is much smaller 
than the transverse trap energy scale, which allows us 
to reduce the problem to one spatial dimension 0, HH . 
The solutions were found numerically using the normal- 
ized gradient flow method [27j, which is able to find a 
state that minimizes the total energy E for given N and 
M, and fulfills the phase matching condition (|12p . The 
conventional imaginary time method is not suitable for 
this problem since the normalization of the wavefunc- 
tion effectively imposes that the chemical potentials p,j 
are equal, thus failing to take into account some of the 
ground states. We note that the direction of the mag- 
netic field with respect to the trap orientation does not 
influence our results, since only the contact interactions 
are taken into account in our model. 

The setup consists of a highly elongated cigar-shaped 
harmonic trap, the longitudinal size of the condensate 
being much larger than other length scales, i.e. the spin 
and magnetic healing lengths £ s and £b. This ensures 
that the contribution of the kinetic energy (e.g. at the 
boundaries between domains) is small compared to other 
energy scales in the system, which is the requirement 
for the application of the Thomas-Fermi approximation. 
Figs.|3[a)-(c) show typical ground state profiles of an an- 
tiferromagnetic condensate in various regimes of param- 
eters. In accordance with the theoretical considerations, 
the domain structure, in general, consists of spatially sep- 
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FIG. 4: Diagrams of phase separation of an antiferromagnetic (a,b) and a ferromagnetic (c,d) condensate in the quasi-lD case. 
The left column corresponds to the untrapped case (spontaneous phase separation), while the right column corresponds to a 
condensate in a harmonic potential. 



arated p$, 2C, and p+ phases, and the arrangement of 
domains is such that the more magnetized the phase, the 
farther it is from the trap center. In the 2C phase, the 
difference between the density of m = +1 and m — — 1 
atoms is constant, in agreement with the condition (1201) . 

As a rule, the 2C+/3+ state (a) is the ground state at 
low magnetic fields, the three-phase state (b) occurs in 
the transitory regime, and the po+P+ state (c) appears 
at stronger magnetic fields. This can be understood by 
noting that the Zeeman energy of the m = component 
is decreased with respect to the m = ±1 components 
when increasing the magnetic field strength. This is the 
reason for the appearance of the po domain which eventu- 
ally replaces the 2C phase completely, in accordance with 
the condition (|2"5|) . Another possibility is the unmagne- 
tized condensate (M = 0), where the whole condensate 
consists of a single po domain (not shown). 

In the case of a ferromagnetic condensate, the ground 
state generally consists of PM and po domains, as de- 
picted in Fig. 03d). The magnetized PM phase is placed 
at the center, again in agreement with the results of 
Sec. [TV] At low magnetic fields the po domains can 
become very small and eventually disappear when their 
size becomes smaller than the spin healing length. In 
this case, the condensate is completely composed of the 
PM phase. In a unmagnetized condensate (M = 0), the 



phase separated PM+po state is preferred at low mag- 
netic fields, SE < 2|c2|n [l6| . otherwise a single po do- 
main is present. 

In Fig. |4l we summarize the above results in a system- 
atic way and juxtapose them with the ones obtained in 
the case of an untrapped condensate [l6[] . It is clear that 
the phase separation has a qualitatively different charac- 
ter in the presence of the trapping potential. Although 
there are similarities, in general the ground states have a 
different structure. In particular, the antiferromagnetic 
condensate can separate into three distinct phases in the 
presence of the trapping potential, while the maximum 
number of phases is two in the untrapped phase. In the 
case of a ferromagnetic condensate, the potential induces 
separation into two phases, while the untrapped conden- 
sate is always homogeneous. 

In Fig. [5] we present an example of a phase separated 
po+p+ ground state of an antiferromagnetic condensate 
in the two dimensional geometry. The domains arrange 
radially, with the unmagnetized phase in the center, in 
agreement with the results of Sec. HVl Atoms of the m = 
— 1 component are present only at the boundary between 
the two domains. We note that in the 2D geometry it 
would be difficult to obtain a clear three phase state, due 
to the limited number of atoms available in experiment. 
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x [|im] 

FIG. 5: The po + p+ state in a quasi-2D harmonic trap. The 
densities of the (a) m = +1, (b) m = and (c) m = — 1 
components are shown, together with the cross section of the 
densities along the x axis (d). The parameters are N = 2.5 x 
10 6 23 Na atoms, m = 0.8, B = 0.15 G, 5£/c 2 n max = 0.15, 
tv xy = 2ir x 10 Hz, tv z = 2ir x 10 3 Hz. 
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FIG. 6: Crossover from the Thomas-Fermi regime to the sin- 
gle mode regime and the disappearance of spin domains. The 
parameters are the same as in Fig. EJb), except that the num- 
ber of atoms was decreased (a) 4x and (b) 16 x , while the trap 
frequency was increased to keep the density constant. Here 
the spin healing length is £ s = 17^m and the magnetic healing 
length is £b = 42^im. 
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FIG. 7: Generation of the ground state profiles by adia- 
batically changing the magnetic field. The magnetic field is 
gradually increased from zero to the final value B = 0.15 G 
during t — 4 s for sodium condensate (a,b), and decreased 
from B = 0.5 G to B = 0.35 G during t = 3 s (c,d) for rubid- 
ium. The frames (a) and (c) show the time dependence of the 
atom density in the m — and m = +1 component, respec- 
tively, and the right column shows the final domain profiles. 
The parameters are as in Fig. [3] except m — 0.6 in (a,b) and 
m — in (c,d). 

transforms gradually towards a single-mode ground state. 
The thickness of the boundaries between domains is de- 
termined by the energy difference between neighboring 
phases. In particular, the boundary between 2C and p+ 
domains is of the order of £ s , while the boundary be- 
tween po and 2C domains is significantly larger, of the 
order of £b- This is clearly visible in Fig.[6jb), where the 
size of the condensate is about two times larger than £ s 
but already smaller than £g. in this case, small p + do- 
mains can still be distinguished, but 2C and po domains 
merged into a single APM phase domain, which is the 
ground state in the SMA limit [H,[l6j]. 



B. Crossover between the Thomas-Fermi and the 
single-mode regime 

The above results, corresponding to the Thomas-Fermi 
(TF) regime, are characteristic for the case when the size 
of the condensate is much larger than the spin healing 
length £ s and the magnetic healing length £g. On the 
other hand, in the opposite limit of a small condensate, 
the ground state structure is well described by the single 
mode approximation (SMA) [lij], where all three mag- 
netic components share the same spatial profile. How- 
ever, typical experimental conditions may correspond to 
an intermediate regime, where condensate size is compa- 
rable to £ s or £g. An example of the crossover between 
TF and SMA regimes is presented in Fig. [5] The param- 
eters correspond to the three-phase state from Fig. [3jb) , 
but the number of atoms was decreased 4x and 16 x 
in (a) and (b), while the trap frequency was increased 
to keep the density constant. The domain structure 



VI. GENERATION OF SPIN DOMAINS 

From the practical point of view, it is important to 
propose a reliable method for the generation of the above 
ground states in experiment. We show that the method of 
adiabatic change of the magnetic field, introduced in [28[, 
can be used for this purpose. In the case of an antifcr- 
romagnetic condensate, we start from the initial ground 
state where all the atoms are in the m = — 1 sublevel. 
Next, part of the atoms is transferred to the m = +1 
sublevel [29| . and the magnetic field is switched off. In 
this way, we obtain a system close to the zero-field ground 
state with an arbitrary magnetization, determined by the 
amount of atoms transferred. Next, we gradually increase 
the magnetic field strength in an adiabatic process ac- 
cording to the formula 



B — \ -Bfinal (26) 
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where t — at the beginning of the switching process, 
^switch is the switching time, and Sg na j is the desired fi- 
nal value of the magnetic field. The form of Eq. (|26|) 
assures that the quadratic Zeeman splitting grows lin- 
early in time. In Figs. [7ta,b) we present an example of 
the generation process simulated by numerical solution 
of the Gross-Pitaevskii equations ([TU]) . Fig. [TJa) shows 
the time dependence of the atom density in the initially 
unoccupied m — component, and Fig. [7th) shows the 
final domain profile. This should be compared with the 
ground state profile in Fig. 02[b), which corresponds to 
the same parameter values. 

In the case of a ferromagnetic condensate, the above 
scenario has to be slightly modified. We begin with 
strong magnetic field, in the regime where the quadratic 
Zeeman energy dominates. The initial state consists of 
atoms in the m = sublevel only, which is the ground 
state in these conditions, see Fig.^Jd). Here, we are lim- 
ited to the states with the total magnetization M = 
only. The magnetic field is then decreased, and after 
crossing the point of phase transition additional PM do- 
mains appear. This scenario and the example of the final 
state are demonstrated in Figs. EJc,d). Since M = 0, the 
curves corresponding to densities of m = — 1 and m = 1 
atoms are overlapping. The PM domain is located in 
the center of the trap, where all three components are 
nonzero. 



VII. CONCLUSIONS 

We investigated spin-1 Bose-Einstein condensates 
trapped in a harmonic potentials under the influence of a 
homogeneous magnetic field. We demonstrated that the 
trapping potential has a strong influence on the structure 
of ground states in these systems, and can make the anti- 
ferromagnetic condensate separate into three, and ferro- 
magnetic condensate into two distinct phases. We stud- 
ied the crossover from the Thomas-Fermi regime to the 
single mode approximation regime, where the conden- 
sate size becomes smaller than the spin healing length 
and the spatial structure of ground states disappears. 
We suggested an experimental method for creation of 
these states by adiabatically changing the magnetic field 
strength. 
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